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Abstract In the dynamics of accretion disks, the pres- 
ence of collective effects associated with the self-gravity 
of the disk is expected to affect not only the momen- 
tum transport, but also the relevant energy balance equa- 
tions, which could differ substantially from the non-self- 
gravitating case. Here we follow the picture that, when 
the disk is sufficiently cold, the stirring due to Jeans- 
related instabilities acts as a source of effective heating. 
The corresponding reformulation of the energy equations 
allows us to: (i) demonstrate how self-regulation can be 
established, so that the stability parameter Q is main- 
tained close to a threshold value, with weak dependence 
on radius; (ii) rediscuss the opacity properties in the self- 
gravitating regime. In particular, we show that, if cool- 
ing is dominated by bremsstrahlung, an optically thin sta- 
tionary accretion solution is thermally stable, even in the 
non-advective case, provided the disk is self-gravitating. 
The details of the cooling function have little effect on 
the structure of such accretion disk, which is in any case 
brought by self-gravity to self-regulate. This condition of 
self-gravitating accretion is expected to be appropriate for 
the outer regions of many disks of astrophysical interest. 
With the reformulation of the energy equations described 
in this paper we have also secured: {Hi) a starting point 
for the study of the emission properties of self-gravitating 
accretion disks; (iv) a tool to analyze the structure of the 
transition region, where the disk becomes self-gravitating. 
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1. Introduction 

A key step in the modeling of accretion disks is the de- 
tailed description of the balance between heating and cool- 
ing terms in the energy transport equations (see Pringie 



1981). In this context, two hypotheses have been generally 



considered. The first states that the only source of heat- 
ing in the disk is viscous dissipation, so that the available 



power is provided by the release of gravitational binding 
energy from the infalling matter. The second hypothesis 
is that cooling is essentially radiative. This latter hypoth- 
esis leads to different disk models, depending on whether 
the disk is optically thick (as in the "standard" model 
of Shakura & Sunyaev 1973) or optically thin (as in the 
model of Shapiro, Lightman & Eardley 1976| ; hereafter 
SLE). One more recent development in the description 
of the energetics of accretion disks has been the recogni- 
tion that radiative cooling may sometimes be insufficient. 
Under these circumstances, most of the energy remains 
stored in the disk rather than being radiated away, lead- 
ing to the so-called advecMon- dominated accretion flows 
(ADAF; Narayan & Yi |1994[) . This situation is described 



by a new equation of energy transport relevant to the ra- 
dial direction, including an advection term (see Sect. |^ 
below) . 

Once the general characteristics of stationary models 
are established, one can proceed to address the problem 
of their stability. Different types of instabilities have been 
considered, leading to important conclusions about the 
possibility of applying the models to concrete systems. In 
one line of research, the stability analysis is carried out 
with the goal of identifying microscopic collective mech- 
anisms able to justify the viscosity a-prescription that 
is commonly introduced, following Shakura & Sunyaev 
( 1973 ), in order to bypass the low level of viscosity pre- 
dicted by classical arguments. Other interesting stabil- 
ity aspects have also been investigated. In particular, the 
thermal-viscous instability related to hydrogen ionization 
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(Meyer & Meyer-Hofmeister 1981) has been invoked to 
explain the variability of dwarf novae. In contrast, the 
optically thin SLE model has been commonly regarded 
as inapplicable, because of the thermal instability that is 
found to occur when optically thin bremsstrahlung is the 
dominant cooling mechanism. To be sure, the thermal sta- 
bility properties of the models depend strongly on the as- 
sumptions about the balance between heating and cooling 
terms (for example, it is well known that, in the ADAF 
context, optically thin solutions exist in the inner disks 
and are thermally stable; see Section || a nd Na rayan & Yi 
1995 ). Indeed, it has been noted (Piran [l97^ ) that small 
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variations in the viscosity law can change the stabihty of 
accretion disk models drastically. 

The self-gravity of the disk can change this picture sig- 
nificantly. It may have an obvious impact o n the radial and 
on the vertical force relations (Paczyhski |l97§| ). Further- 
more, it has long been recognized that non-axisymmetric 
instabilities associated with the self-gravity of the disk can 
affect the properties of angular momentum transport and 
hence the relevant viscosity la w (e.g ., see Lin & Pringle 
199C, Laughlin & Bodcnhcimer 1994), although it still re- 



mains to be clarified whether such modifications can be 
reconciled with the standard a-prescription. Recently we 



have argued (Bertin 1997, Bertin & Lodato 1999) that 
self-gravity should be properly incorporated also in the 
energy balance equations, which should thus take a more 
general form than usually considered. In these first in- 
vestigations we actually replaced the energy equations by 
a self-regulation prescription demanding that the axisym- 
metric stability parameter Q is maintained close to a value 
Q basically dictated by the condition of marginal Jeans 
stability. The physical justification for the existence of an 
efficient dynamical thermostat of this kind is based on the 
presence of suitable competing heating and cooling mech- 
anisms, but was not made explicit in terms of a set of 
consistent equations for the energy transport. This sim- 
plified approach is adequate for the construction of dy- 
namical models, and has already led to some unexpected 
results (in particular, the natural occurrence of flat rota- 
tion curves). Yet, it deliberately avoids the issue of what 
is the predicted emission of the accretion disk, which has 
had such an important role in the confrontation of models 
with observations in this general research area. We antici- 
pate that the disk self-gravity should be significant only in 
the cold outer disk. In this respect, the identification of a 
consistent set of energy equations would also be desirable 
in order to set up a proper matching of solutions across 
the relevant transition region. 

In this paper, we include the effect of Jeans-related in- 
stabilities in the energy balance equations by means of a 
term of effective heating. Such heating term is modeled 
as a function steeply rising when Q falls below a thresh- 
old value Q and practically vanishing when Q rises above 
Q. We then show how this more general energy equation 
leads to the self-gravitating, self-regulated accretion disk 
solution, characterized by a value of Q very close to Q 
and only weakly dependent on radius. This stationary so- 
lution, although non-advective and optically thin, turns 
out to be thermally stable. 

The paper is organized as follows. In Section ^ we 
briefly summarize the arguments that lead to the balance 
of heating and cooling terms in current theories of accre- 
tion disks. In Section ^ we propose a way to incorporate 
the effects of self-gravity in the energy equations by means 
of an additional heating term. In Section ^ we then provide 
and discuss a complete set of equations for the construc- 
tion of the full models of self-gravitating accretion disks. In 



Section ^ we outline the structure of the stationary models 
and study their thermal stability. The conclusions of the 
paper are given in Section ^, where we briefly comment 
on possible astrophysical applications. 



2. Energy balance in non-self-gravitating disks 

A generally accepted approach to the construction of ac- 
cretion disk models starts from the assumption that the 
disk is geometrically thin (so that vertically integrated 
and vertically averaged quantities can be usefully intro- 
duced) and adopts a fluid model, specified by the continu- 
ity equation, the radial component of the Euler equation, 
the conservation of angular momentum in the presence 
of viscous torques, a viscosity prescription (such as the 
a-prescription), an expression for the vertical thickness 
of the disk (to be derived from the solution of the ver- 
tical hydrostatic equilibrium equations), and the relevant 
energy equations. In the following, we will refer to verti- 
cally integrated or averaged quantities, such as the surface 
density a, as "surface" quantities, while we will refer to 
other quantities, such as the volume density p, as "vol- 
ume" quantities. 

The energy equations are generally based on two hy- 
potheses, of purely viscous heating (with the addition, in 
some cases, of irradiation from the inner disk or from the 
central object) and of purely radiative cooling. The first 
hypothesis is thus contained in an expression for the heat- 
ing rate in terms of "surface quantities" : 



Q+ = lyairQ'f = D{r), 



(1) 



where v is the viscosity coefficient and il = ri(r) is the 
differential rotation (usually taken to be Keplerian). 

The second hypothesis, that cooling is entirely radia- 
tive, implies that the cooling rate per unit surface Q~ can 
be calculated from the "volume" properties of the disk 
and the assumed radiation mechanisms: 



Q- ^R. 



(2) 



The energy equations are then split into an equation 
for the radial transport and an equation for the verti- 
cal transport. In practice, in systems with given physical 
properties, one should carefully check the validity of this 
separation, which is natural only in the limit of thin disks. 
For the horizontal transport the equation is the so-called 
"advection" equation: 



auTQ— = — Q' 
ar 



(3) 



where u is the radial velocity, Tq is the midplane temper- 
ature, and s is the entropy per unit mass. This equation 
states that, locally, the unbalance between energy gained 
from heating and energy lost from cooling is stored in 
the disk and transported radially along with the infalling 
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matter. Entropy in Eq. (|^) is calculated based on the first 
principle of thermodynamics: 



ms„=dg = c„dT + pd(l/p), 



(4) 



where p is the volume density, p is the pressure, is the 
specific heat at constant volume and s„ is the "volume" 
entropy. After averaging over the vertical direction one 
defines the averaged "surface" entropy so that: 



J^O-T- — — 

ar 7 



1 dc2 



dpo 

1 dr po dr ' 



(5) 



where c is the thermal speed on the equatorial plane, po is 
the central volume density, so that a = 2pQh {h is the ge- 
ometrical thickness) and 7 is the specific heat ratio. Note 
that in Eq. one neglects the effects of horizontal ra- 
diative energy transport (but see Eq. (3) in Narayan & 
Popham 1993). For the vertical transport one uses the 
equations of radiative transport, which, integrated over 
the z-coordinate, give rise to Eq. (||). 

A final remark is in order. The above equations are 
a set of differential equations, and hence they should be 
supplemented by a set of appropriate boundary condi- 



tions. However, for standard disks (see Sect. ^ below), in 
which one neglects advection terms, the system reduces to 
a set of algebraic equations and no boundary conditions 
are thus required. In contrast, ADAF solutions are based 
on the full differential problem. Even in this case, however, 
the self-similar solution found by Narayan & Yi (1994) 



(see Sect. 2.2) avoids the need of specifying the boundary 
conditions. Other ADAF models take such self-similar so- 
lution as the approximate solution to be matched at large 



radii. In this context, it has been noted (Yuan 199£) that 
the boundary conditions are, indeed, important in deter- 
mining the structure of the accretion disk and should be 
handled carefully. 

2.1. Standard disk 

The so-called standard disk solutions are those for which, 
in Eq. (^, the dominant terms are and Q^, so that 
the advection equation reduces to: 

Q+«Q-, (6) 

If we express this relation in view of the physical hypothe- 
ses mentioned previously, then the energy balance equa- 
tion for standard disks is: 



D = R. 



(7) 



As R describes the emission of radiation from the disk, 
Eq. ^ shows that the emission is related to in a simple 
way, largely independent of the detailed construction of 
the model. On the other hand, Eq. (^, involving all the 
surface quantities (ct, Tq,...), is the closure relation needed 
to calculate the structure of the accretion model. 



2.2. ADAF models 

In some cases, the output of radiation may be insufficient 
so that the balance between and Q~ in Eq. (^ is 
not possible. In this more general case one cannot use Eq. 
(^ (and hence Eq. (0)). This situation has often been 
described in terms of a convenient parameter, called /, 
defined as: 



/ 



so that Eq. (||) can be rewritten as: 



ar 



(8) 



(9) 



For specified radiation processes, it is clear that Eq. 
(||) and Eq. ^ together can close the model, which is 
characterized (for given values of the viscosity parameter 
a and of the mass and angular momentum accretion rates 
M, j) by a profile / = /(r) determined by Eq. (H). In 
reality, it has been noted that, at a given radius r, by 
considering the relation between / and M (a and J being 
fixed), there are in general three branches of solutions (at 
least for sufficiently small values of M): (i) one branch 
corresponding to / <C 1, so that (see Eq. (^)) R « D, 
is similar to the solution of standard disks, (ii) one new 
branch, called ADAF, with / « 1, is characterized by the 
balance between viscous heating and the advection term 
(see Eq. (H)); for this branch there is little emission of 
radiation, as i? = (1 — f)D ^ D. Finally, a third branch 
(iii), for intermediate values of /, is thermally unstable 
and corresponds to the SLE solution. 



Treating / as a given parameter, Narayan & Yi ( 1994 ) 
have shown that models closed only by Eq. (^ (that is, 
not taking into account the radiation properties) admit a 
simple, self-similar solution. This solution can be seen as a 
first approximation to the ADAF solution on a wide (but, 
in any case, limited) radial range. 

3. Influence of self-gravity on the energy balance 

From the solution of the non-self-gravitating disk equa- 
tions, it is interesting to check the behavior of the param- 
eter Q, which, as is well known, regulates the stability of 
the disk against axisymmetric Jeans instability. For a fluid 
one-component disk the parameter Q is defined as: 



CK 



(10) 



where k is the epicyclic frequency (equal to f2 in the Kep- 
lerian case). Self-gravity becomes important at large radii. 
In fact, Q{r) is expected to be a decreasing function of ra- 
dius. For example, for the optically thick disk dominated 
by the Kramers opacity (the outer region of the Shakura & 
Sunyaev 1973| disk) it is proportional to r^^/^. When the 
standard model is applied to concrete cases (for example, 
the conditions typical of an AGN) Q is found to plunge 
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well below unity in regions of interest (Kumar 1999). This 
behavior is untenable. It is only a reminder that the equa- 
tions considered initially to construct the accretion disk 
model must be reformulated so as to incorporate the im- 
portant role of the disk self-gravity (e.g., see Kolykhalov 



& Sunyaev 1980) 



As a quick, constructive way to address the problem 
of self- gravitating accretion disks, we have recently argued 



(Bertin 1997) that the energy equations, when self-gravity 



is dominant, should be replaced by the physically based 
prescription that Q is self-regulated to a value Q close 
to unity. This is all we need if we are interested in cal- 
culating the main dynamical characteristics of such self- 
regulated accretion disks. On the other hand, a detailed 
examination of the physical justification at the basis of 
the self-regulatio n pro cess (see Sect. 3 of the article by 
Bertin & Lodato 1999 ) readily shows that self-regulation 
results from heating and cooling mechanisms that go be- 
yond those considered in Sect. ^. In the self- regulation 
picture, widely considered in the study of spiral galaxies, 
the Jeans instability has the effect of inducing a fast heat- 
ing of the disk (on the dynamical time-scale), at a rate 
that depends strongly on Q. Such heating thus tends to 
bring Q towards values of stability. On the other hand, if 
the system is Jeans stable to begin with, efficient cooling 
mechanisms (that, in some cases, may be just the radia- 
tive cooling of Eq. (||)) can cool the disk down to lower 
values of Q. In the following Sections we will reformulate 
the energy balance equations starting from the arguments 
just presented. Such reformulation of the energy equations 
will allow us to: {i) demonstrate how self-regulation can 
be established; (n) rediscuss the opacity properties in the 
self-gravitating regime. In addition, it will give us {iii) a 
starting point for the study of the emission properties of 
self-gravitating accretion disks; {iv) a tool to analyze the 
structure of the transition region, where the disk becomes 
self-gravitating. A study of the last two points is post- 
poned to a separate investigation. 

3.1. Additional heating 

We suggest that in the expression of the heating rate a 
term Hj should be added (so that (5+ = D+Hj), strongly 
dependent on the value of Q, which we express in terms 
of surface quantities as: 



Hj = g{Q)ac^n. 



(11) 



Here g{Q) is a function that rises quickly when Q drops 
below Q (for example, g{Q) ~ (Q/Q)", with n ^ 1). 
The form of this heating term is argued on dimensional 
grounds, given the fact that the collective effects associ- 
ated with self-gravity are expected to heat the disk on the 
dynamical time-scale. The general structure of Hj is anal- 
ogous to that of the heating term introduced in the context 
of spiral structure in galaxies to demonstrate the time- 
dependent process of self-regulation (see Bertin & Romeo 



1988 ). The argument is also analogous to that used by 
Lin & Pringle (1987) to incorporate the collective effects 
associated with self-gravity into an effective viscosity. In 



fact, we could write (if we take h 
Hj = {ajch)a{rn')^, 
with the definition: 



'-/nGa): 



g{Q) 



(12) 



(13) 



With this device, we see that the term we are introduc- 
ing mimics the viscous heating for a system with modified 
viscosity a ^ a + aj. In practice, we think that the two 
issues, of a modified viscosity in the presence of spiral 
instabilities and of a heating term associated with Jeans- 
related instabilities, should be kept separate. Note that 
in the context of spiral galaxies the emphasis is reversed, 
and little attention is paid to the problem of viscosity, 
although it is well known that spiral instabilities eventu- 
ally favor mass accretion toward the center. In any case, 
it should be stressed that the present approach remains 
heuristic. The issue of energy dissipation requires further 
clarification (see also Balbus & Papaloizou 1999| ). 



3.2. Additional cooling 

By referring to the context of galactic disks, we find an- 
other important element that is generally overlooked in 
the discussion of the energy balance equations for accre- 
tion disks. 

Basically this is the distinction between thermal speed 
and effective thermal speed. In practice, the effective 
thermal speed entering the dynamical equations may be 
higher than the thermal speed entering the radiation equa- 
tions. The interstellar medium of galaxy disks is a cloudy 
medium, for which the thermal speed associated with the 
temperature of the clouds is much smaller than the ve- 
locity dispersion of the cold gas clouds. It is the latter 
that sets the appropriate effective thermal speed relevant 
to the dynamics of the disk. When a fluid model of the 
interstellar medium is adopted, in studies of disk dynam- 
ics, the numbers that are used for the equivalent acoustic 
speed are then those applicable to the turbulent speed 
characteristic of the cloudy medium. Therefore, one im- 
portant cooling mechanism is that associated with the in- 
elastic collisions between gas clouds. Mechanisms of this 
kind may operate also in accretion disks, but they are 
usually ignored (but see Quataert & Chiang 20001 ). 



3.3. Modified energy balance equation 

Because of these considerations we see that the energy bal- 
ance equation for a non-advective disk, should be written 
as: 



D + Hj>R. 



(14) 
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In the inner disk Q is likely to be much larger than 
unity: if we drop the caveat about the distinction between 
thermal speed and effective thermal speed, we see that 
such equation would then reduce back to that of a stan- 
dard disk. 

4. Basic equations for steady-state self-gravitating 
accretion disks close to the transition region 

In the following we will study the problem of self- 
gravitating accretion disks in a region not far from 
the transition region from non-self-gravitating to self- 
gravitating disk (see Fig. 

Equations 

The analysis is thus simplified in two respects. 

(j) On the one hand, the radial balance of forces for a 
cool, slowly accreting disk will be described by the Kep- 
lerian relation: 



(15) 



where is the mass of the central object. In Eq. 
( |l5| ) we have neglected the contribution of the disk self- 
gravity. This approximation is valid for r <C = 
2{GMi,){Q/Af(GM /2a)-2/3, where a is the viscosity pa- 
rameter (see Eq. ( pO| ) below). For parameters typical of 
an AGN, previous studies (Bardou, Heyvearts & Duschl 
1998| ) show that the transition from non-self-gravitating 



to self-gravitating disk occurs at a radius, which we call 
rg, that is much smaller than r^, so that, at least for a 
certain parameter range, the approximation in Eq. ( p^ ) is 
applicable to the inner parts of the disk. The effects asso- 
ciated with the integral relation between the gravitational 
field generated by the disk d$cr/dr and a are described in 
the paper by Bertin & Lodato ( 1999] ). 

(zz) On the other hand, the energy balance equation is 
taken to be (for the definitions of D and Hj see Eqs. (0) 
and (ITT 



(16) 



D + Hj^R 



with R specified by the relevant radiation mechanism. 
In the following we will use the expression for the cool- 
ing rate per unit volume appropriate for optically thin 
bremsstrahlung (but see the additional discussion of Sect. 
5j): 



2'KkT 2*6^ 
3me ShnieC^ 



gn 



1.410-2^r^/^gn2, 



(17) 



where n is the particle density, given by n = p/jirup (fi 
is the mean molecular weight), and ^ « 1.2 is a gaunt 
factor averaged over frequencies (see Rybicki & Lightman 
1979| ), and where the last expression is evaluated in CGS 



units {T measured in K). In Eq. (|17| ) c is the speed of 
light and should not be confused with the effective thermal 



speed defined above; we have assumed that the atomic 
weight is unity, and that rii = Ug. In the spirit of the 
thin disk approximation we will define R = qbr{z = 0)h. 
The optically thin description breaks down at small radii, 
below a scalelength r^, where the disk is optically thick 
(see discussion in Sect. |5|). 

The remaining equations are the equation of continu- 
ity: 



M 



-2Trrua, 



the conservation of angular momentum: 
J = Mr^Vt + 2iTvar^ — , 



(18) 



(19) 



and the standard a-prescription (Shakura & Sunyaev 
1973) for viscosity: 

V = ach, (20) 
where h is the vertical thickness of the disk, defined by: 



t:G(j 4Q2 



v/l + 8Q2/7r-l 



(21) 



The last equation is an interpolation formula for the "ex- 
act" thickness calculated including both the self-gravity of 
the disk and the gravitational force of the central object in 
the vertical hydrostatic equi libriu m (see Eq. (A9) in the 
article by Bertin & Lodato 1999 and discussion below), 
in the limit of Keplerian rotation curve. Note that in the 
limit of non-self-gravitating disk (Q S> 1) this equation 
leads to the commonly used result h = ^ 7r/2(c/f2). In the 
following we will take J = 0; the effects of a non-zero net 
angular momentum flux modify the disk structure only in 
the innermost parts of the disk. 

The level of viscosity present, as measured by the pa- 
rameter a, should depend also on the development of spi- 
ral instabilities, and, therefore, on the parameter Q. Since 
for the self-gravitating part of the disk Q will turn out 
to be self-regulated, we prefer to simplify the analysis by 
leaving a as an assigned parameter. Our reference value 
(often considered for AGN-type accretion disks) will be 
a = 0.05. 

4.2. Characterization of the transition region 

The above set of equations is naturally reduced to dimen- 
sionless form by the introduction of the length-scale: 



TT _ 4 / GM 
3 ^ hc^ 9a \ 3a 



-1/3 



{GAQ_ 
2-kG 



1/2 



2/3 



.(22) 



In practice, this length-scale roughly corresponds to the 
out ermost radius for which our solution exists (see Sect 
|3.l| below) . In the following we will thus refer to the di- 
mcnsionless radius f — r /rj^. 
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deviations from Keplerian rotation 





"standard" disk 
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transition from 
standard disk 
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Figure 1. Sketch of the disk structure and of the different 
zones with the corresponding scale lengths. In practice, the 
present paper addresses the third zone, between rr and Vs- 



Based on the scale lengths introduced so far, and on 
general results that will be better justified in Sect. ||, we 
can characterize the structure of the accretion disk, in 
the region we are focusing on, by the following four zones 
(see Fig. P: (?) An inner standard disk with the value of 
Q decreasing rapidly outwards; (ii) A transition zone, in 
which the heating term due to self-gravity becomes impor- 
tant and the optical thickness of the disk decreases until 
(Hi) the optically thin regime is reached; (iv) At radii of 
the order of Vs self-gravity affects also the rotation of the 
disk and at even larger radii the simple set of equations 
adopted in this Section breaks down and possibly other 
cooling mechanisms come into play. 

4-.3. Additional remarks on the role of self- gravity in the 
region considered 

The term "self-gravitating" is often used with different 
meanings, because the quantity irGa may be judged to be 
important in relation to different quantities: (a) The cu- 
mulative effect of irGa may make the gravitational force in 
the plane different from the Keplerian field generated by 
the central point mass; in this respect, to check whether 
the disk should be considered to be self-gravitating, we 
should compare the location r of interest with the scale 
Ts of the problem. In the limit r 3> the gravitational 
field of the central point mass may be neglected, but at 
r w the disk is already self-gravitating, (b) The trans- 
port properties in the disk may be affected significantly by 
density waves and other Jeans-related phenomena; here 
the importance of the disk self-gravity is quantified by 
comparing irGa to the acceleration ck and thus it is quan- 
tified by Q. Note that, in this respect, for a disk to be self- 
gravitating it is sufhcient that Q be close to unity; in other 
words, the self-gravitating limit is not set by the condition 
Q <C 1, as might be thought at first, (c) The vertical scale 
of the disk may be significantly changed with respect to 
that of a non-gravitating disk; as shown in a separate ar- 



ticle (Bertin & Lodato |l99g| ; see Fig. Al and Eq. (A9) 
in the Appendix), this depends both on the local value of 
Q and on the overall mass distribution in the disk. Fur- 
thermore, the expression of the disk thickness given in Eq. 



(A9) of Bertin & Lodato (1999) is valid uniformly and can 
be considered as an "exact" formula; the self-gravitating 
expression h = (? /ttGct is applicable, for example, when 
Q 1 and the gravitational field in the disk plane is 
dominated by the contribution of a disk characterized by 
u ^ \ jr. 

To a large extent, the three items above are indepen- 
dent aspects of the role of self-gravity, even though we 
may argue that at very small radii any accretion disk is 
likely to be non-self-gravitating from all the three points 
of view, while a cold, radially extended accretion disk at 
large radii is bound to be self-gravitating with respect to 
all aspects (a), (b), and (c). 

In this paper we are less interested in the specific im- 
pact of item (a) (investigated earlier by us) and wish to 
focus on item (b), which had not been properly explored 
in relation to the energy equations. Fortunately, effects re- 
lated to Jeans instability, which we model by means of the 
heating term described by Eq. (|ll|), become important at 
a location r « tq, which, for models applicable to AGN 
disks, occurs well far in with respect to the location where 
item (a) comes into play, because rq ^ rg. As to the disk 
thickness, we have just recalled that the expression given 
in Eq. (A9) of Bertin & Lodato (1999) can be considered 
as an "exact" formula; thus its Keplerian limit, as adopted 
in Eq. (|2l|), is legitimate when item (a) is unimportant. 
Therefore, in the transition region we can safely follow the 
set of equations described in subsection 4.1. 



5. Thermal stability of the self-regulated solution 

Optically thin, non-advective disks have been generally 
discarded because of the thermal instability that results 
from an increasing of the heating rate relative to the cool- 
ing rate in response to an increase of the temperature. This 
occurs because they are found to violate the condition for 
thermal stability: 



>o, 



(23) 



where the subscript a indicates that the derivative should 
be taken at constant surface density (see, for example, 
Pringle [19811 ). 

Since our discussion of Sect. ^ modifies the basic in- 
gredients of the energy balance equations, we should now 
check what are the properties of the self-gravitating solu- 
tions in relation to the thermal stability condition. Given 
the fact that the additional heating term Hj given in Eq. 
(|ll[) is a function of Q and that such parameter is pro- 
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portional to the thermal speed c, we can reformulate the 
thermal stability condition as: 



dQ 



Log 



> 



(24) 



icy 
o 




Figure 2. Relation between the ratio of heating to cooling 
rates and the stability parameter Q at three different radii, 
i.e., from right to left f = 10"'', 10~^, 1. The stable branch 
is the one with positive slope. The intersection between the 
vertical line and the stable branch is the stable self-regulated 
solution. 



5.1. The thermally stable self- gravitating branch 

In Fig. H we show the relation between /Q^ and Q 
at fixed radius, for different values of the radius. Equilib- 
rium corresponds to Log{Q~ / Q^) = (the vertical line 
in Fig. ^) and the solution is thermally stable if the slope 
of the curve is positive. As is easily recognized, there are 
two branches of solution: one, unstable, is characterized 
by Q 3> 1, for which the dominant heating mechanism is 
viscous dissipation as in non-self-gravitating disks (this so- 
lution can be regarded as analogous to the SLE solution). 
The other solution is self-gravitating (it has Q « 1) and it 
is, indeed, thermally stable. Note that such stable equilib- 
rium solution is always characterized by a value of Q close 
to unity. The equilibrium Q increases slightly with radius 
(although this dependence is weak); independently of the 
physical parameters of the disk, the maximum value of Q 
for the stable solution is w 1.3Q. There is a maximum ra- 
dius for which this kind of solution is possible, located near 
tr. The relative location of the two scale lengths rR and 



defines the importance of the deviation from Keplerian 
rotation in our model. For parameters typical of AGN con- 
figurations {AU = 10^Mq,M = IM^/yr « 610^^g/s) 
we estimate rn/rg ~ 510^. Hence, before the scale is 
reached the rotation curve is expected to depart from the 
simple Keplerian curve, so that the quantitative results 
of the present simplified work are only approximate when 
f > 2 10~^ and may even be qualitatively incorrect when 
f is of order of unity. The lack of solutions at large radii 
if > 1) is due to the fact that radiative cooling is in- 
sufficient and unable to balance the heating rate. In this 
context, on e should keep in mind that, as we indicated 
in Sect. 3.2, other kinds of effective cooling could be effi- 
cient and play a role in the outer regions, such as inelastic 
collisions between clouds, which have been neglected here. 

5.2. Consistency of the optically thin description 

We have also calculated the optical depth of the disk, 
to check that the model of optically thin disk is consis- 
tent. We express the " effective" optical depth r* (Ry- 
bicki & Lightman 1979 ) for the case kt ^ Hbr (where 
K.T and Kjjr are the opacity due to Thomson scattering 
and bremsstrahlung, respectively), as: 

(25) 



T* = \/nrTT, 

where Tbr = KbrO' and tt = ktO'. 




0.001 



Figure 3. Optical thickness of the disk as a function of radius 
for M = IMQi/yr. In the innermost regions the disk is optically 
thick. 



In Fig. ^ we show the optical depth as a function of 
radius for the thermally stable equilibrium solution. The 
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Figure 4. Ratio of Thomson opacity to free-free opacity, 
showing that Thomson scattering always dominates over 
bremsstrahlung. 



disk is optically thin in the outer regions and the opti- 
cal depth increases inwards. In the innermost regions the 
disk eventually becomes optically thick. This occurs at a 
dimensionless radius tt independent of the central mass 
and very weakly dependent on the mass accretion rate (ap- 
proximately fr M~^/^^); for M = IMQ/yr the transi- 
tion to the optically thick region occurs at an inner radius 
fr = tt/tr « 410~^. In Fig. ^ we show the profile of 
the ratio kt / Hbr to check that Thomson opacity always 
dominates over free-free absorption. 

In the transition region from optically thick to opti- 
cally thin emission, for a proper matching, a more detailed 
description would be desired. 

5.3. Robustness with respect to changes of the radiative 
processes 

In contrast with other cases for which optically thin emis- 
sion is considered (i.e., the inner regions of disks with very 
low mass accretion rates), in our case the temperature of 
the disk (for reasonable values of M) is expected to be low, 
so that pure bremsstrahlung fails to be the dominant ra- 
diative process and line emission and absorption become 
important. In view of this remark, we have checked the 
consequences of a different cooling function, using a func- 
tional interpolation (Sarazin & White 1987) of numerical 
data (Raymond, Cox & Smith 1976), appropriate for lower 
temperatures. 

In this test of robustness, the value of the accretion rate 
M, through the temperature of the outer disk defined as 



kTout = {iJ.mp/2){GM /Sa)"^^^ , appears explicitly in the 
problem. This is the result of the fact that the various 
line emissions become efficient at specified temperatures. 
In practice, only minor differences are found for the self- 
gravitating branch (even under extreme conditions; see 
Fig. H), with respect to the previous analysis, because the 
cooling function influences only the equilibrium value of 
Q but this is, in any case, self-regulated to values close to 
unity (on the thermally stable equilibrium solution). 




bremsstrahlung 

T„„, = 7*10"K 
T„^^=3.5*10'K 
T =7*10''K 



-0.5 

2 4 6 8 

Log(QVQ') 

Figure 5. Stability curves for different choices of the cooling 
function at radius f — 10~^. The details of the cooling mech- 
anism do not change the main characteristics of on the stable 
self-gravitating branch of solutions. The three outer tempera- 
tures correspond to M ~ 90, 30, 3MQ/yr, respectively. Here 
the value a = 0.05 has been assumed. 



6. Discussion and conclusions 

Self-gravitating accretion disks may be present in many 
astrophysical systems, ranging from protostellar disks to 
AGN and could also be of interest for the general galactic 
context. 

In protostellar disks it is likely that, especially at the 
beginning of the accretion processes, during which the 
mass of the central object is small, self-gravity plays an im- 
portant role. The influence of self-gravity has been some- 
times invoked to explain the flat infrared spectrum of 



some T Tau stars (Adams, Lada & Shu 1987, Duschl, 
Strittmatter & Biermann 2000). However, current esti- 



mates for the mass accretion rates in these systems (typ- 
ically M ~ 10~^ MQ/yr) so far have discouraged inter- 
pretations of the shape of the spectrum in terms of self- 
gravitating, optically-thick disks. 
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On the other hand, observational evidence suggests 
that, in some cases, the currently accepted models may be 
inadequate to describe the accretion flow of some AGN. 
For example, in the case of the Seyfert galaxy NGC1068, 
deviations from Keplerian rotation have been observed; 
when the standard Shakura & Sunyaev ( 1973| ) model is 
applied to this case, the values of the stability parameter 
Q in the outer disk turn out to be exceedingly low (Kumar 
1999). In fact, the energy balance in standard disks, for 



which viscous dissipation is exactly balanced by radiative 
heating, leads to a Q profile rapidly decreasing outwards. 

In this paper, we have focused on the problem of the 
energy budget in self-gravitating accretion disks, with the 
ultimate goal of providing a framework to describe the 
emission and the spectrum associated with such disks. We 
have adopted the view that the effects of self-gravity on 
the energy balance are not only related to viscous dis- 
sipation, but that there are also processes, such as the 
self-regulation of Jeans instability, that should give a sig- 
nificant contribution to the heating of the disk, at least in 
its outer parts. 

By referring to conditions applicable to some AGN, we 
have then shown that a thermally stable, non-advective, 
optically thin equilibrium solution is possible in the region 
dominated by self-gravity. In contrast with other optically 



thin solutions (Shapiro, Lightman & Eardley 1976), our 
solution is thermally stable as a result of the strong de- 
pendence of the heating term related to self-gravity on 
the parameter Q. We have also found an unstable solu- 
tion, which is not self-gravitating (Q 3> 1), and we argue 
that it is a modification of the SLE solution. 

The value of Q on the thermally stable solution is al- 
ways close to unity, and it is slightly dependent on radius. 
The temperature in the outer disk is constant and de- 
pends on the mass accretion rate M . The structure of the 
disk, as long as the energy balance is dominated by self- 
gravity, turns out to depend only weakly on the detailed 
radiation processes involved. This has been shown by tak- 
ing into account two different radiation processes: pure 
hremsstrahlung and the case of a cooling function domi- 
nated by line absorption and emission (Sarazin & White 



1987| ), which is best suited for the low temperature con- 
ditions expected in the outer disk. The dependence on 
the radiation processes is only weak because of the self- 
regulation mechanism: whatever the radiation process, the 
disk readjusts, so as to keep Q close to marginal stability. 

For simplicity, in this paper we have focused on condi- 
tions relevant to the region of the disk where the transi- 
tion to the self-gravitating regime occurs. Further out, one 
should include the presence of additional cooling mecha- 
nisms, not considered here, and especially the effects re- 
lated to the deviations from a Keplerian rotation curve. 
The latter effects have b een w idely discussed in a previous 
paper (Bertin & Lodato 199£ ). 
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